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Context. The problem of the existence of intermediate-mass black holes (IMBHs) at the centre of globular clusters is a hot and 
controversial topic in current astrophysical research with important implications in stellar and galaxy formation. 
Aims. In this paper, we aim at giving further support to the presence of an IMBH in u> Centauri and at providing an independent 
estimate of its mass. 

Methods. We employed a self-consistent spherical model with anisotropic velocity distribution. It consists in a generalisation of the 
King model by inclunding the Bahcall-Wolf distribution function in the IMBH vicinity. 

Results. By the parametric fitting of the model to recent HST/ACS data for the surface brightness profile, we found an IMBH to cluster 
total mass ratio of M./M = 5.8^' x 10~ 3 . It is also found that the model yields a fit of the line-of-sight velocity dispersion profile 
that is better without mass segregation than in the segregated case. This confirms the current thought of a non-relaxed status for this 
peculiar cluster. The best fit model to the kinematic data leads, moreover, to a cluster total mass estimate of M = (3.1 ± 0.3) x 10 6 M Q , 
thus giving an IMBH mass in the range 1.3 x 10 4 < M. < 2.3 x 10 4 M Q (at l<x confidence level). A slight degree of radial velocity 
anisotropy in the outer region (r ^ 12') is required to match the outer surface brightness profile. 



Key words, black hole physics - stellar dynamics - methods: analytical 
- globular clusters: individual: ll> Centauri (NGC 5139) 



methods: numerical - galaxies: kinematics and dynamics 



1. Introduction 

Intermediate-mass black holes (IMBH) still belong to the class 
of 'exotic' objects in the current astrophysical belief. With 
masses between M. ~ 100-10 4 M o , they would represent the 
minor mass counterpart of super-massive black holes - whose 
existence is established with much more robustness - but still 
more massive than stellar black holes. One of the places where 
they should more likely be located is among the densest stellar 
environments in the Universe, i.e. at the globular clusters (GCs) 
centre. However, so far, the most direct observable signature of 
their existence, namely the emission in the radio and X-ray bands 
(mainly from Bondi-Hoyle accretion o f intracluster gas), is not 
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compatible with the i nferred presence of a ~ 6 x 10 3 M IMBH 
(lLanzoni et al ., 2007). On the other hand, no detectable X-ray 
sources have been found at the centre of mass of NGC 2808, 
leading IServillat et all d2008l) to state that M. £ 290 M in this 
cluster. 

In fact, in most cases only upper limits for IMBHs 
masses can be deduced fro m radio observations (see, e.g. 
iMaccarone & Servillall I2008L for NGC 2808 and for a gen- 
eral discussion). These surp risingly low upper limits lead 
IMaccarone & Servillatl (|2008) to cast doubts on the fact that the 
scaling relation M. - cr, where cr is the central velocity dis- 
persion of the host stellar system (with mass M and luminos- 
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yet really clear and conclusive (see, e . g..lLiu & Di Stefanol2008t ltv L )> is the same M - ~ o" 4 S law that has been clearly noted 



2009t andlMiller & Colberdl2004[ for a general review). 

To date, only the GC Gl (in M31) exhibits a detected source, 
see n in both radio (with a n 8.4 GHz power of 2 x 10 15 W Hz -1 , 
see lUlvestad et al.Ll2007h an d X-ray (with a 2 x 10 36 erg s -1 lu- 
minosity at 0.2-10 keV, see iPoolev & Rappaporl 120061) bands. 
The observed fluxes, as well as their ratio, are compatible with 
the cl aimed presence of a ~ 2 x 1O 4 M IMBH (Geb hardtet all 
120051) - altho ugh other kinds of sources cannot be completely 
ruled out (e.g. lKong et all 120091) . Another extra-galactic hyper- 
luminous X-ray source (5 - 100 x 10 40 erg s" 1 at 0.3 - 10 keV) 
is located in the SO-a galaxy ES0243-49 and its features suggest 
an IMBH emission. Recently, an unresolved optical counterpart 
with brightness comparable t o that of a massiv e GC has been 
identified around this source (So ria et al ., 2009), though higher 
resolution observations are needed. 

In our Galaxy, the central region of NGC 6388 hosts an unre- 
solved s et of X-ray sources, with a total luminosity of 2.7 x 10 33 
erg s^ 1 dNucita et all [2008), implying an accretion efficiency 



for super-massive black hol es in galaxies dFerrarese & MerrittT, 
2000; Gebhar dt et al l l2000l) . It should be emphasised, however, 
that the upper limits on IMBH masses drawn from X-ray or ra- 
dio observations strongly depend on the assumption that the gas 
distribution around the compact object is uniform (isotropic ac- 
cretion). It is clear that, if this distribution had any amount of 
dumpiness, those limits could be largely underestimated. 

Nevertheless, the question of the validity of the extrapola- 
tion of the M. - cr scaling relation to IMBHs is still open and 
deserves to be discussed briefly here. In general, this relation 
can be understoo d as a consequence o f the fundamental scal- 
ing law M. oc M dMagorrian et all 119981). In galaxies, this scal- 
ing law and the two relations M ~ L 5 ^ 4 dFaber et all 1 19871) and 
L~ cr 4 dFaber & JacksonLll976h . lead just to M. ~ cr 5 . In glob- 
ular clusters, on the other hand, the obser ved trends are M ~ L 
and L ~ cr 5/3 dMevlan & HeggieL Il997l) . which in fact yield 
M. ~ cr . This implies a generally lowe r mass ratio between the 
IMBH and the host cluster, as noted by IMaccarone & Servillatl 
(2008]). A shallow M. - cr relation, namely M. ~ cr 12 , was al- 
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ready reported in Miocchi (2007, hereafter M07) based on para- 
metric IMBH mass estimates (see below). On the other hand, we 
must mention a recent study of this specific topic, in which the 
low-mass extrapolation of the galactic M. ~ cr 4 8 relation seems 
to fit a sample of 5 reported IMBHs in GCs ( Saf onova & Shastrl 
l20loh . 

In view of all this, the study of possible IMBH fingerprints 
on either star-count or surface brightness (SB) profiles is a detec- 
tion route that still deserves to be pursued, especially when kine- 
matic observations close to the IMBH gravitational influence re- 
gion are available. In this respect, a spherical and self-consistent 
model of GCs with a central IMBH at rest was presented in 
IM07L both with equal mass stars, i.e. the single-mass (SM) case, 
and with a multimass (MM) stellar spectrum inclunding mass 
segr egation. This model is an extension of King-Michie mod- 
els dMichid 119631: iMichie & Bodenheimerl 119631; iKind, 119661) 
that is obtained by including the Bahcall-Wolf stellar distribu- 
tion function within the IMBH gravitational influence region. 
The latter was shown to solve the Fokker-Planck equation in 
the vicinity of a ce ntral IMBH that formed long before a clus- 
ter re laxation time (Bah call & Wohj 1 19761: iBinnev & Tremaina 
1 1987b . and its validity was subsequently confirmed by accurate 
numerical simulations (Freitag & Benz. 2002; Baumgardt et all 
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The typical SB profile that comes out of the model has, for 
any reasonable IMBH mass, the appearance of a normal low- 
or medium-concentration cluster (c ;$ 2) and shows a steep 
cusp only in the very inner region (typically within a tenth of 
the core radius) delimited by the 'cusp radius', r cu . In fact, out- 
side r cu , a shallow power-law behaviour - with a logarithmic 
slope s $ 0.25 - is the most easily observable fingerprint in 
the otherwise flat core profile. Interestingly, this confirmed the 
finding of other authors who - using a completely different ap- 
proach (accurate A^-body simulations) - also claim that IMBHs 
most likely reside in non-core-collapse d clusters showing just a 
weak rise of the S B in the core region (Baumg ardt et al.l l2005; 
iTrenti et aDl2007l) . Recently, high-resolution Montecarlo simu- 
lations have provid ed another inde pendent confirmation of these 
structural features dUmbreit et al.| [2009). On the other hand, ac- 
cording to other A^-body experiments, it is claimed that post- 
core-collapsed GCs also exhibit a King-like profile, but with a 
s ~ 0.4 - 0.7 steep core beha viour dTrenti et all 120101) ; it must 
be emphasised, however, that lM07l models yield core profiles 
that are always significantly flatter and unable to fit behaviours 
with such a high s. 

The shape of the SB profile given by the M07 model de- 
pends on 2 dimensionless parameter^ For the purposes of this 
study, we use the IMBH to cluster mass ratio, M./M, and the 
dimensionless gravitational potential at the edge of the IMBH 
dynamical influence region, Wbh- The latter replaces the usual 
King model's central dimensionless potential Wo, with the aim 
of avoiding the si ngula rity at the cluster centre in the presence 
of the IMBH (see lM07l for further details). 

In[M07 it was shown that lower and upper limits of M. exist 
as a function of c and s. This relationship was then applied to in- 
vestigate the presence of IM BHs in the set of GCs, who se SB 
was accurately measured in iNovola & Gebhardl d2006l) using 
HST/WFPC2 archive images. Among the six candidate clusters 
found, NGC 6388 and M54 have subsequently been the objects 
of further and more detailed studies (through parametric fitting 



of star-count profiles) that suggest the presence of an I MBH with 
mass ~ 6 x 1 3 Mp in the forme r dLanzoni et alll2007h and ~ 10 4 
in the latter dlbata et all 120091 in this case kinematic data were 
also exploited). On the other hand, the massive cluster u> Cen 
was not chec ked as a possible can didate, because it was not in- 
cluded in the lNovola & Gebhardtl d2006l) sample, and moreover, 
small slopes in t he core region could not be revealed in published 
SB profiles (e.g. lMevlariri987l;lFerraro et alll2006l) . 

Nonetheless, a recent and accurate determinati on of the a> 
Cen c entre and the use of HST/ACS images led INovola et all 
(2008) to detect a steeper profile in the core region of this pe- 
culiar cluster, thus suggesting the influence of an IMBH. By fit- 
ting the high inner peak of the line-of-sight velocity dispersion 
(LOSVD) found from Gemini GMOS/IFU integrated light spec- 
troscopy (= 23 + 2 km s" 1 at an average radius ~ 1".9) with non- 
para metric and orb it-based models with uniform mass-to-light 
ratio. iNovola et al J estimate a ~ 4 x 10 4 M Q object residing at the 
cluster centre. However , by solving the spherical an d aniso tropic 
Jeans equation on the lAnderson & van der Mare ll (20101) pro- 
j ected density and kinematic data. Ivan der Marel & Anderson I 
d2010h find that the presence of an IMBH is possible only if 
M./M ^ 4.3 x 10~ 3 , wh ich corresponds to about half the mass 
predicted by Novo la et all 

In this paper we would like to provide further evidence on 
the presence of the IMBH and to give another independent es- 
timate of its mass, by means of a parametric fitting of both the 
SB and the LOSVD profiles using the M07 model. The results 
from the best fit of the SB profile are described in Sect. [2] while 
those coming from the LOSVD fitting are presented in Sect. [3] 
Concluding remarks are reported in Sect. [4] 



2. The fit of the surface brightness 

To study the a> Cen SB profile , we considered the HST/ACS 
measurements recently made by INovola et al.l d2008l) inside 40" 
from the c luster centre, w hile for outer radii we took the obser- 
vations by Mevlanl dl987L their Table 1). 

As is evident from the too low concentration of the dotted 
profile in Fig. [TJ (bottom panel), we notice that the SM isotropic 
model is unable to fit the outermost part (r ^ 13') of the SB 
profile, where, however, the background contamination should 
be negligible, for it was shown to only be relevant for r k, 33' 
dLeon et alll2000l) . Thus, the discrepancy from the prediction of 
this model should be due to the intrinsic dynamical state of the 
cluster outskirts. In fact, we find that a good fit of the whole 
profile can be achieved by including either a certain degree of 
radial velocity anisotropy or an MM stellar population with mass 
segregation (keeping isotropic velocities). 

Nevertheless, the fit of the inner SB profile with a SM 
isotropic model permitted us to determine the best fit value for 
M./M regardless of velocity anisotropy, because the presence of 
the latter only influences the ou ter region (as happens i n normal 
King- Michie models, see, e.g., iGunn & Griffinl [l 9791 : iMiocchil 
2006). Thus, a grid of SM isotropic models have been generated 
by sampling the form parameters Wbh and M./M. As the model 
profiles are expressed in dimensionless units, they have to be 
scaled in both the radial and the SB dimension. Thus, for each 
model of the grid we found the best fit values for two suitable 
scale parameters, namely the "visual" core radius^ r c and the 



1 The model can also include velocity anisotropy in the GC outskirts. 
In this case the outer SB profile shape depends on the anisotropy radius, 
too. 



2 Here r c is defined as the radius at which the SB drops to half 
its value at r cu . In good approximation, r c coincides with the loca- 
tion of the "turn-off" of the profile (also called 'break radius' in 
Novola & Gebhardt 2Q0|; see lM07l for more details). 



P. Miocchi: A mass estimate of an intermediate-mass black hole in a> Centauri 



3 



1 \ ' -I 
i 




log(r[arcsec]) 

Fig. 1. Bottom panel: surface brightness profile of a> Cen and 
the SM isotropic (dotted line) and anisotropic (solid line; with 
r a = 4.5r c ) best fit. In the case of the isotropic model, only data 
with log(r) < 2.9 have been considered for the best fit search. 
The MM isotropic model yields a best fit profile that is indis- 
tinguishable from that given in the SM a nisotropic ca s e. For 
log(r) < 1.6, the HST/ACS observations bv lNovola etall (120081) 
are u sed (filled squa res), while for log(r) > 1 .6 data are taken 
from Meylan (1987, open squares). The top panel shows an en- 
larged view of the central region. The core radius is plotted with 
a dashed line. 

normalisation value of the SB, restricting the fit to data points 
with r < 13'. 2. Since the SB measurements uncertainties are 
known, we minimised the x 2 sum weighted by the width of the 
error bars. 

The calculated^ 2 values are reported in Fig. [2] from which 
we deduce that W BH = 5.25 and M./M = 5.8^ x 1(T 3 , with a 
level of confidence (LOC) of 68.3%. The best fit isotropic model 
gives r c = 156", which conf irms both the more rec ent observa- 
tions b y lFerraro et"aT1 ((2006) and the value listed in lTrager et alj 
(119951) . 

Once Wbh and M./M has been determined by fitting the 
inner SB profile, we fit the entire data set by including radial 
velocity anisotropy in the st ellar system o utskirts (confirming 
what was already noted by Mevlan 1987), namely outside an 
anisotropy radius r a ; see, e.g., lMiocchiT ?2006) for a description 
of how anisotropic velocities can be efficiently implemented in 
King-Michie models. Thus, a "sub-grid" of anisotropic models 
is generated by sampling r a in the range [2, 10] x r c . The result- 
ing x 2 behaviour, this time evaluated over all SB data, is plot- 
ted in Fig. [3] and leads to the estimates r a = (4.5 + 0.1) x r c = 
12' + 0'.2 with an LOC of 68.3%. From this figure, it can also 
be noted how the anisotropic model when 'pushed' towards the 
isotr opic case (r a » r c ) giv es unacceptable fits (huge x 1 val- 
ues). Ivan de Ven et alJ (|2006) find that the velocity distribution 
in a> Cen is nearly isotropic inside ~ 10', in agreement to our 
best fit value for r a . On the other hand, these authors reported 
the presence of a slight tangential anisotropy in the cluster out- 
skirts which is, however, below the uncertainty in the velocity 
dispersion measurements (see their figure 8). 




Fig. 2. Contours of x 2 as a function of Wbh and of the M./M 
ratio, plotted for the grid of SM isotropic models (the filled dots, 
with size proportional to^ 2 ) fitted to data points with r < 13'. 2. 
The central cross marks the minimum^ 2 model location (Wbh = 
5.25 and M./M = 5.8 x 10~ 3 ) and the isocontours for Ay 2 = 2.3, 
6.17, 1 1.8 correspond to confidence regions of 68.3 (thick line), 
95.4, and 99.7%, respectively. 




4.2 4.4 4.6 4.8 

r / r c 



Fig. 3. x 2 behaviour of the best fit SM model in the isotropic 
case (the cross in Fig. |2]i when anisotropy is introduced, as a 
function of the anisotropy radius r a . In this case, all the available 
SB data are taken into account. The x 2 values (for Ax 2 = 1, 
9) corresponding to an LOC of 68.3 and 99.7% are also plotted 
(dotted lines). 



The SM anisotropic best fit model reported in Fig. [TJ yields 
a tidal radius r t = 41 '.3 and a concentration parameter, c = 
log(r t /r c ) - 1-2, su bstantially lower than the 1.6 value quoted 
in the lHarrisJ (|1996[) catalogue, but in good agreement w ith re- 
cent results dFerraro et all 2006; va n de Ven et all 2006). 

As we said at the beginning of this section, the entire ob- 
served SB profile can al so be fi t ted fa irly well by an MM 
isotropic model (see also iMevlanl 1 19871) including the IMBH 
and stars in the mass range 0.4 — 1.2M Q distributed following 
the Salpeter mass function, with central energy equipartition. 
In this case, the best fit profile is for M./M = 2.8 x 10~ 3 
(and Wbh = 7.5), and it practically overlaps with that of the 
anisotropic case (Fig. Q]). Nonetheless, we discarded this MM 
model because it underestimates the LOSVD in the central re- 
gion, as we see in detail in next section. 
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3. Velocity dispersion profile 

To provide an estimate of M., we have to quantify the cluster 
total mass M first. This can be done by exploiting the most re- 
cent kinematic observations of a> Cen. For this purpose, along 
with the two innermo st points taken from the Gemini GMOS- 
IFU measuremen ts in | Novola et"a"D ( 20081) . we use the LOSVD 
data employed by van d e Ven et al. see references therein 

for the discussion of the various data sources). These are based 
on various independent sets of measurements, which in most of 
the radial annuli include values taken in different apertures (see 
Fig.|4]i. A radial error bar is plotted for the innermost point to in- 
dicate the width of the 5" x5" GMOS-IFU field of view that was 
centred on the cluste r nucleus to obtain the integrated spectrum 
dNovola et all 120081) . 

Once the SB profile has been fitted, the form of the LOSVD 
profile is univocally given by the model and cannot be adapted 
to the observed behaviour. The best fit can be found by adjusting 
only the velocity scale factor (corresponding to a vertical shift- 
ing of the profile). In turn, this factor depends on the adopted 
cluster distance and total mass M. We chose to a-priori fix the 
distance to 4.8 kpc (as from lvan de Ven et alll2006h . and then to 
find the M value that gives the best fit to the LOSVD observa- 
tions. Shown in Fig.|4]are two LOSVD best fit profiles: the one 
given by the SM anisotropic model and the one produced by the 
MM isotropic one. 

It is evident (see also Fig. [5J that the SM case yields a better 
fit to LOSVD data, having P(x 2 > x\ x ) = 84%, compared with 
the MM model that gives P(x 2 > Xa) = 47% mainly because, 
in the inner region (log r < 1 .8), it exhibits too low an LOSVD. 
This is naturally expected from the mass distribution in the MM 
case being dominated by the lighter (and fainter) stars, which 
are much less concentrated than the giants. Thus, the velocity 
dispersion of the giants starts to decrease at larger radii, con- 
sequently the best fit tends to give a lower inner LOSVD in the 
attempt to fit the outer data. Interestingly, that the SM model bet- 
ter represents the dynamical situation of this cluster suggests that 
mass segregation has not been efficient in a> Cen. In this sense, 
it confirms the current thought that this cluster is not completely 
relaxed by coll isions, because of its relatively l ong relaxation 
time (see, e.g.. lMevlanlll987t iMevlan et al.lll995l) . Various au- 
thors, indeed, have found indications of a uniform mass-to -light 
ratio (see, e.g JMerritt et all 1 19971: Ivan de Ven et al.Ll2006l) . 

Considering the relatively large error bars of LOSVD mea- 
surements in crowded regions, from Fig. [4] we note that the 
model profile predicts that the cen tral velocity c u sp is a pparently 
more centrally concentrated than Nov ola et alJ {2008) observa- 
tions suggest (it starts to be evident only for r % 1"), though one 
has to consider the "visual effect" of the logarithmic scale in r. In 
fact, a relatively large residual (~ 4 km s _1 ) still remains for the 
innermost LOSVD data point, although it s radi al error bar inter- 
sects the model profile (at r - 0".6). If the lM07l model represents 
the real cluster dynamical state well, this could indicate the in- 
fluence of some statistical bias affecting this bin or too large an 
average radius chosen for it. In this respect, it is also worth not- 
ing that recent and accurate proper motion measurements reveal 
no significant velocity cusp at the central region of this clus- 
ter, though this study relies on a different dataset and kinematic 
centre location and, moreo ver, the authors do not observe an y 
appreciable cusp in density dAnderson & van der Marell20loh . 

In iNovola et alJ (120081) the innermost bins are fitted quite 
well (apart from the measurement at log(r) 1.7, see their 
Fig. 4). However, it must be noticed that in these authors' model 
the IMBH mass best fit value depends almost completely on the 



few innermost LOSVD data points, while it has practically no ef- 
fects for r t 30" and plays no role at all on the SB profile. In our 
case, on the contrary, the behaviour of the LOSVD given by the 
model is strongly dependent on the best fit parameters of the SB 
profile. If our model were forced to fit the entire LOSVD well, 
then the required ~ 4 x 1O 4 M IMBH would produce a much 
steeper SB core behaviour (along with too low a concentration), 
which would be completely different from the observed one. The 
disadvantage of our parametric approach is that it is "less gen- 
eral", because it is constrained by the theoretical hypothesis ly- 
ing behind the assumption of that particular distribution function 
in phase-space. 

The predicted cluster total mass is M = (3.1 ± 0.3) x 10 6 M 
with an LOC of 68.3% (Pig. [5j. It is in marginal agreement with 
the dy namical estimate of (2.5+O.3)xlO 6 M in lvan de Ven et all 
(120061) - thou gh it would agr ee well at 2<x level - while much 
lower than fhe lMevlanl (1 19871) 3.9 x 1O 6 M value. This author 
used a King-Michie MM model in which an approx imated en- 
ergy equipartition was imposed (see Miocchi, 2006, for a dis- 
cussion of this approximation). This, together with the assumed 
presence of very low-mass stars (down to 0.13 M ), can explain 
the higher M estimate. As far as the M/L ratio is concerned, if 
one assumes a total V-band lumin osity in the "p r udential" range 
L v = (1.0 ± 2) x 10 6 L (e.g. ISeitzeri 119831: iMevlanl Il987t 
Carraro & LiaL 120001) . one gets M/ L v = 3.1 ±0.9, a value com- 
patible with the accurate 2.5 + 0. 1 Ivan de Ven et al] J2006) esti- 
mate. Of course, as the model predicts no mass segregation, the 
mass-to-light ratio turns out to be uniform. 

The estimate made in Sect. |2]of the ratio M./M, combined 
with the uncertainty on the cluster total mass, yields an IMBH 
mass in the range 1.3 x 10 4 < M. < 2.3 x 1O 4 M (with a 
68.3% LOC ), which span s abou t one third to a half the mass 
predicted by INovola et al.l d200 8 |). but is marginally comp atible 
with the $ 1.3 X 10 4 M^ lvan der Marel & An derson (2010) esti- 
mate. However, it has to be kept in mind that these two estimates 
rely on different cluster centre. Finally, it is worth noting how 
our estimate range, though still incompatible, gets closer to the 
~ 2300 M upper limit as constrained by the oj Cen radio con- 
tinuum emission dMaccarone & Servillatl 12008). 



4. Conclusions 

In this paper we have presented a parametric fit of the sur- 
face brightness (SB) profile of a> Cen (N GC 5139), made up 
of HST/A CS data i n the central region dNovola et all [2008) 
and of the IMevlanl (1 19871) normalised profile in the outskirts. 
The fit was done by using a self-consistent (spherical and non- 
rotating) model that includes a central intermediate-mass black 
hole (IMBH). The whole SB profile (from r ~ 1".6 out to r ~ 
42') can be well-fitted by the model both with the single-mass 
stellar distribution - assuming a radially anisotropic velocity dis- 
tribution outside 12' - and with a multimass model with isotropic 
velocity. The comparison of the generated LOSVD with the 
kinematic observations recently enriched at the very central 
regio n by Gemini GMOS-IFU measurements dNovola et all 
2008), however, allows this degeneracy to be resolved in favour 
of the single-mass case. In fact, the multimass model yields 
too low an LOSVD in the central region. This suggests that a> 
Cen is presently in a non m ass-segregated state, as already ar- 
gued by various a uthors (e.g. lMevlan1ll987llMevlan et all ll 995: 
Ivan de Ven et all 120061) . It is also worth noting that recent N- 
body studies show that the presence of an IM BH in a cluster can 



body studies show that the presence of anlfy lt 
suppress mass segregation dGill et all 12008). 
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Fig. 4. Best fit LOS VD radial profile in the SM anisotropic (solid 
line) and MM isotropic case (dashed line; in this case, it is de- 
termined by weighting the contributions of all the stellar com- 
ponents according to their luminosity). The best fit total cluster 
mass is 3.1 and 4.4 x 10 6 M o , respectively . Filled squares are 
the m easurements fr om Novo la et alJ (|2008), open squares come 
from lvan de Ven et al.l d2006l) . 




M [10 6 M s ] 

Fig. 5. x 1 behaviour of the LOSVD fit as a function of the to- 
tal cluster mass M in the SM anisotropic (solid line) and MM 
isotropic (dashed) case. The dotted lines indicate the 68.3 and 
99.7% LOC. 

From this parametric study we deduce the 68.3% confidence 
intervals M = (3.1 ± 0.3) x 10 6 M for the cluster total mass and 
M./M = 5.8^ 2 x 10 3 for the mass ratio, leading to an estimate 
for the IMBH mass in the range 13,000 < M. < 23,000M o . 
This value is from abo ut one third to a half the ~ 40, 000 M Q 
mass predicted bv lNovola et alJ ([2008), though it is compatible 
with th e 18, 000 M n upper limit provided b y the dynamical anal- 
ysis in Ivan der Marel & Anderson Note, however, that 
these two published estimates are based on different cluster cen- 
tre. 
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